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PREFACE 
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Reclamation,  Agreement  No.  9-AA-60-00830.  Appreciation  is  expressed  to  Darryl  Calkins 
and  Rae  Melloh  for  helpful  discussions. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes. 
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Ice  Effects  on  Hydraulics  and  Fish  Habitat 


GEORGE  D.  ASHTON 


INTRODUCTION 

\; 

In  the  temperate  zones  of  the  world  the  formation  of  ice  in  streams  and  rivers  may  have 
significant  effects  on  their  water  depths  and  velocities,  even  when  the  stream  discharge  does 
not  change.  Velocity  and  water  depth  are  among  the  attributes  that  constitute  the  habitat  of 
the  aquatic  species  that  reside  in  the  streams  and  rivers.  Thus,  to  evaluate  the  habitat  it  is 
important  to  be  able  to  characterize  the  effects  of  ice  on  velocity  and  depth. 

In  this  report  I  summarize  the  effects  of  river  ice  on  hydraulic  behavior,  with  examples 
meant  to  provide  guidance  in  evaluating  what  habitats  may  be  expected  in  winter  in  com¬ 
parison  with  non-ice  conditions  of  the  same  stream.  While  river  ice  formations  may  be  com¬ 
plex  at  the  small  scale  of  tens  of  meters  or  less,  on  a  larger  scale  they  are  often  sufficiently 
uniform  to  enable  reasonable  comparisons  of  depth  and  velocities  with  and  without  ice 
present.  The  emphasis  here  is  on  shallow  rivers,  such  as  the  Platte  River  in  Nebraska.  The 
Platte  is  subject  to  many  competing  demands  for  its  flow,  among  which  are  the  demands  for 
habitat.  A  concern  for  the  required  winter  flows  was  the  motivation  for  this  study.  ^  <*  / 


THEORETICAL  BACKGROUND 
Steady  flows  without  ice 

For  a  stream  of  uniform  slope  S,  the  mean  velocity  V  (feet  per  second)  is  related  to  the 
hydraulic  radius  R  (ft)  (defined  as  the  cross-sectional  area  divided  by  the  wetted  perimeter) 
by  the  Manning  equation: 

V  =  141.  R2/i  s  1/2  (1) 

n 

where  n  is  an  empirical  (Manning’s)  coefficient  associated  with  the  roughness  of  the  bed 
surface.  By  multiplying  by  the  cross-sectional  area  A,  this  may  be  put  in  the  form  of  a  dis¬ 
charge  relation: 

Q  =  AV  =  142-  AR  2/3  S  1/2  (2) 

n 

where  Q  is  the  discharge.  Most  streams  are  much  wider  than  they  are  deep,  so  the  hydraulic 
radius  is  nearly  equal  to  the  depth  {R  =  D),  and  the  cross-sectional  area  A  =  BD,  where  B  is 
the  width  and  D  is  the  depth.  For  non-ice-covered  flow,  then, 

Qo=  Av  =  1AS_  BDo  D2/3  Som 

nb 


(3) 


Table  1.  Values  of  Manning's  n.  (After  Rouse  1950.) 


0.016-0.017  Smooth  natural  earth  channels,  free  from  growths, 
with  straight  alignment. 

0.020  Smooth  natural  earth,  free  from  growths,  little  cur¬ 

vature.  Very  large  canals  in  gixxl  condition. 

0.022  Average,  well-conslircied.  moderate-si/ed  earth 

canals  in  good  condition. 

0.025  Very  small  earth  canals  or  ditches  in  goixl  condition. 

or  larger  canals  with  some  growth  on  hanks  or  scat¬ 
tered  cobbles  in  bed. 

0.030  Canals  with  considerable  aquatic  grow  th.  Rock  cuts. 

based  on  average  actual  section.  Natural  streams  with* 
good  alignment,  fairly  constant  section.  Large 
lloodwav  channels,  well  maintained. 

0.035  Canals  halfchoked  with  moss  growth.  Cleared  but  not 

continuously  maintained  flood  ways. 

0.040-0.050  Mountain  streams  in  clean  loose  cobbles.  Rivers  with 

variable  section  and  some  vegetation  growing  in 
banks.  Canals  with  very  heavy  aquatic  growths. 

0.050-0. 1 50  Natural  slreamsof  varying  roughnessand  alignment. 

The  highest  values  for  extremely  bad  alignment,  deep 
pools  and  vegetation,  for  floodways  with  heavy  stand 
of  timber  and  underbrush. 


where  the  subscript  o  denotes  the  open  water  case  and  »()is  the  roughness  of  the  bottom.  The 
open  water  case  is  shown  in  Figure  1  a.  The  value  of  the  Manning  roughness  factor  is  generally 
determined  empirically  from  measurements  of  R,  S  and  T  in  the  field;  however,  there  is  suf¬ 
ficient  experience  that  a  good  estimate  can  be  made  from  general  descriptions  of  the  boundary 
surface.  Table  1  gives  typical  values  for  various  surfaces  ranging  from  concrete  to  earth  with 
weeds.  Typical  values  are  on  the  order  of 0.02-0.04  for  most  natural  streams  and  higher  for 
mountain  streams.  Note  also  that  n  has  dimensions  of  ft 1/6  to  make  eq  1  dimensionally 
homogeneous,  and  the  coefficient  1.49  has  dimensions  of  R  (g)  where  ,t>  is  gravitational 
acceleration  (m  s-*). 


Steady  flow  with  ice 

A  similar  analysis  of  steady  flow  in  a  stream  of  uniform  slope  with  ice  may  also  be  done 
(Fig.  1  b).  Since  the  wetted  perimeter  is  now  doubled.  R  =  D./2  and  the  ice-covered  equivalent 
of  eq  3  from  eq  2  becomes 


Qi=  L49  BD.  (Oif35.i/2  (4) 

«c  \  2  I 

where  the  applicable  roughness  value  n  has  a  subscript  c  representing  a  "composite" 
roughness,  that  is,  one  that  depends  on  both  the  roughness  of  the  bed  and  the  roughness  of 
the  underside  of  the  ice.  Much  work  has  been  done  to  determine  the  appropriate  value  of  the 
composite  roughness  as  a  function  of  the  roughness  of  the  bed  and  the  ice  undersurface.  One 
widely  accepted  formula  is  the  Sabaneev  expression; 


« c  =  »b 


i  +  K/»hP 


12/3 


(5) 


Figure  1 .  Schematic  diagram  of  flow. 


where  11  is  the  roughness  coefficient  of  the  ice  undersurf  ace.  For  a  smooth  ice  undersurface 
such  as  often  forms  when  the  ice  thickens  in  place  during  the  winter,  /;.  can  be  as  small  as  0.01 , 
while  for  an  ice  cover  formed  from  fragments  of  broken  ice.  /rcan  be  as  large  as  0.06  or  more. 
The  value  of  n.  for  ice  covers  formed  from  broken  ice  seems  to  depend  on  the  thickness  of 
the  accumulation  making  up  the  ice  cover.  Table  2 


presents  typical  data  (Nezhikovskiy  1964). 

Equation  4  does  not  directly  include  the  ice  cover 
thickness,  only  the  depth  of  flow  beneath  the  ice  cover. 
The  water  level,  of  course,  is  higher  since  the  floating 
ice  cover  displaces  water.  The  height  of  the  water 
(above  the  bed)  is  denoted  by  H  to  distinguish  it  from 
the  flow  depth  under  the  ice  cover.  Since  most  ice 
covers  in  medium-sized  and  large  streams  are  floating, 
the  displacement  effect  is  associated  with  the  specific 
gravity  of  ice  relative  to  water  and 

H  =  0.916  a  +  D.  (6) 

where  t.  is  the  ice  cover  thickness,  and  0.9 1 6  is  the  ratio 

I 

of  the  density  of  ice  to  the  density  of  water. 


Table  2.  Manning’s  n  for  ice  ac¬ 
cumulations  at  freeze-up.  (Af¬ 
ter  Nezhikovskiy  1964.) 


Thickness  of 
accumulation 
(m) 

Manning's 

n 

Loose 

slush 

Dense 

shah 

Ice 

floes 

0.1 

_ 

_ 

0.015 

0.3 

0.010 

0.013 

0.04 

0.5 

0.01 

0.02 

0.05 

0.7 

0.02 

0.03 

0.06 

1.0 

0.03 

0.04 

0.07 

1.5 

0.04 

0.06 

0.08 

2.0 

0.04 

0.07 

0.09 

3.0 

0.05 

0.08 

0.10 

5.0 

0.06 

0.09 

— 

COMPARISON  OF  OPEN- WATER 
AND  ICE-COVERED  CASES 


Equations  3  and  4  provide  a  basis  for  comparing  the  depths  and  velocities  for  the  open- 
water  case  and  the  ice-covered  case,  I  will  assume  that  in  both  cases  the  discharge  is  the  same, 
since  it  generally  depends  on  things  upstream  of  the  reach  of  interest.  I  first  assume  the 


3 


( 


uniform  flow  case  for  which  the  slope  of  the  water  surface  is  the  same  as  the  slope  of  the 
stream  bottom  and  hence  the  same  for  both  open  and  ice-covered  conditions.  From  eq  3 


D0  = 


«b  Q 


13/5 


1.49  55  1/2  J 


(7) 


and  from  eq  4 


D{  = 


r.  »c  g  2* ; 

1 1.49  55  l/2. 


3/5 


(8) 


The  ratio  of  depth  under  the  ice  to  open-water  depth  is.  under  the  conditions  I  have  assumed 
(same  Q,  6  and  S), 


2l 

D0 


.-,2/3 


L«b 


3/5 


(9) 


Neglecting  the  effect  of  n ^/^(discussed  later),  eq  9  shows  that  for  a  constant  slope  the  depth 
under  the  ice  increases  32%  over  the  open-water  case.  Since  the  velocity  is  inversely 
proportional  to  the  depth,  if  the  discharge  is  the  same  ( Q  =  VQD0-  V.Dj),  then  the  velocity 
under  the  ice  is  reduced  by  24%  relative  to  the  open-water  case.  This  effect  is  that  due  to 
doubling  the  wetted  perimeter. 


COMPOSITE  ROUGHNESS  EFFECTS 


I  now  examine  the  effects  of  composite  roughness.  Incorporating  the  Sabaneev  formula 
(eq  5)  into  eq  9  results  in 


£l 

D0 


1.32 


1  +  (n,/nb) 


3/2 


2/5 


(10) 


Figure  2  shows  a  plot  of  D./D^as  a  function  of  The  range  plotted  for  «  /«  is  somewhat 

large  but  met  with  occasionally.  For  example,  a  very  smooth  ice  cover  with  n.  =  0.01  and 
/ib=  0.03  causes  the  depth  to  increase  by  7%,  while  a  rough  ice  cover  with  n.  =  0.06  and 
n  =  0.03  causes  the  depth  to  increase  by  70%. 

Similarly,  for  the  same  discharge,  the  velocity  will  change  from  the  open- watercase  to  the 
ice-covered  case  according  to  the  ratio 


=  D<l  =  0.76 

V'o  D, 


1  +  («,/«b)3/2 


2/5 


(11) 


In  summary,  if  the  slope  remains  constant,  the  depth  will  increase  and  the  velocity  will  de¬ 
crease  as  long  as  n./nb  >  0. 

In  the  above  analysis  I  have  assumed  that  the  slope  is  constant  for  the  two  cases.  For 
shallower  streams  this  is  more  likely  the  case.  However,  particularly  for  deeper  streams,  it 
is  also  possible  for  the  energy  slope  to  increase  locally  with  little  change  in  depth.  I  now 
examine  the  change  in  velocity  when  the  depth  is  forced  to  remain  constant.  From  eq  3 


5 


O  — 


Q  "b 

1.49  5D0V3 


and  from  eq  4 


(12) 


4 


Figure  2.  Effect  of  the  ratio  of  ice  cover  roughness  to 
bed  roughness  on  the  flow  depths  under  the  ice  cover 
for  a  river  of  constant  slope. 


Figure  3.  Effect  of  the  ratio  of  ice  cover 
roughness  to  bed  roughness  on  the  slope  of 
a  river  if  the  depth  is  held  nearly  constant. 


[  1 .49  BDff-  | 

Hence,  for  the  same  discharge,  width  and  depth. 


(13) 


(14) 


After  the  Sabaneev  formula  for  n^  is  applied.  S,/So  is  plotted  in  Figure  3. 

Now,  what  does  this  slope  change  mean?  Suppose  a  stream  or  river  has  a  slope  of 0.000 1 . 
Increasing  the  slope  by  a  factor  of  2.5  (the  case  of  n.  =  «b)  changes  the  slope  to  0.00025.  This 
means  the  slope  can  change  over  a  mile  with  a  depth  change  of  0.5  feet  at  the  end  of  the  reach. 
For  deeper  rivers  and  canals  with  controls  at  the  ends  of  the  reach,  this  is  the  usual  effect,  while 
for  shallower  streams  the  depth  adjusts. 

In  actual  cases,  of  course,  both  the  slope  and  the  depth  adjust,  with  the  relative  amounts 
of  each  determined  by  the  end  conditions  (the  ‘‘boundary  conditions”).  While  somewhat 
oversimplified,  a  generalization  is  that  deep,  shorter  channels  adjust  to  ice  by  changing  slope, 
while  shallower,  longer  channels  adjust  to  ice  by  changing  depth. 


EFFECT  OF  ICE  ON  STAGE 

The  analysis  presented  above  describes  the  change  in  the  depth  of  flow,  i.e..  the  distance 
between  the  undersurface  of  the  ice  and  the  river  bottom.  However,  what  is  usually  observed 
is  the  stage  of  the  river,  i.e..  the  elevation  of  the  top  of  the  water,  or  how  high  the  water  rises 
in  a  hole  cut  in  the  ice  cover.  The  stage  with  an  ice  cover  is  the  depth  of  the  flow  plus  the 


5 


displacement  due  to  the  ice  cover.  This  displacement  depends  on  the  density  of  the  ice.  The 
stage  at  any  point  is  thus 


H  =  D  +  fj  (15) 

Pw 

where  t  =  ice  thickness 
1 

p.  =  density  of  the  ice 
p  =  density  of  the  water. 

For  most  ice  covers  p/pw=  0.916.  Thus,  the  thicker  the  ice.  the  higher  the  stage,  even  with 
no  change  in  the  depth  of  flow  beneath  the  ice  cover. 

OTHER  EFFECTS 

The  analysis  presented  above  is  somewhat  idealistic  because  it  assumes  a  simple 
uniformly  thick  ice  cover  and  straight-sided  channel  cross  sections,  and  it  does  not  ac¬ 
commodate  the  boundary  conditions.  There  are  many  other  effects  that  influence  the  flow 
beneath  an  ice  cover.  Some  of  these  are  discussed  here. 

Blockage  of  the  flow  cross  section  by  ice  deposits 

In  many  natural  rivers,  ice  is  produced  as  frazil,  which  consists  of  small  crystals  of  ice  that 
form  in  the  flow.  Frazil  is  often  seen  in  the  form  of  floating  slush.  This  ice  is  easily  transported 
by  the  flow  and  is  deposited  in  slower  portions  of  the  cross  section.  These  deposits  may  extend 
over  the  entire  depth,  leaving  flow  passages  that  occupy  only  part  of  the  cross  section  beneath 
the  top,  solidly  frozen  ice  cover.  Experience  suggests  that  there  is  a  critical  velocity  below 
which  such  frazil  is  deposited  and  above  which  any  frazil  deposits  are  eroded.  The  actual 
value  of  this  threshold  velocity  is  uncertain,  but  it  seems  to  be  between  2  and  3  feet  per  second 
(fps).  This  suggests  that  for  streams  steep  enough  to  have  velocities  over  about  2  fps.  a  flow 
area  will  be  maintained  to  carry  the  discharge  at  that  velocity  and  the  unneeded  portions  of 
the  cross  section  will  be  filled  with  frazil  deposits.  The  2-fps  value  is  also  approximately  the 
threshold  below  which  a  surface  ice  cover  can  form  and  above  which  arriving  ice  is 
transported  beneath  the  ice  cover,  thus  preventing  upstream  progression  of  the  ice  cover.  If 
the  initial  open-water  velocity  is  above  2-3  fps,  the  ice  production  will  accumulate  and  cause 
a  phenomenon  termed  “staging,"  i.e.,  the  accumulation  raises  the  water  level  and  slows  the 
velocity,  and  the  ice  cover  then  progressively  accumulates  upstream.  The  result  is  a  cross 
section  largely  blocked  by  frazil  deposits  with  channels  through  or  beneath  with  a  character¬ 
istic  mean  velocity  of  2-3  fps. 

Effect  of  partial  ice  coverage 

Often  the  ice  cover  does  not  cover  the  entire  channel.  The  effect  is  to  raise  the  stage  and 
shift  the  distribution  of  flow  such  that  more  occurs  in  the  open  surface  portions  of  the  cross 
section  and  less  in  the  ice-covered  portions.  The  stage  change  is  gradual  as  the  coverage  goes 
from  zero  to  complete.  This  case  is  analyzed  in  Appendix  A. 

Ice  effects  on  multiple-channel  flow  distribution 

In  rivers  with  multiple  channels,  one  of  the  channels  may  be  ice  covered  while  the  other 
remains  ice  free.  The  effect  is  the  same  as  with  partial  ice  coverage  of  a  single  channel.  More 
flow  is  carried  by  the  ice-free  channel  and  less  by  the  ice-covered  channel  than  if  both  were 
ice  free.  The  effect  on  the  smaller  channel  can  be  significant.  This  situation  is  analyzed  in 
Appendix  B  for  the  case  of  two  channels.  The  presence  of  hydraulic  controls  other  than  ice 
may  constrain  the  ice  effects. 
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APPLICATION  TO  THE  PLATTE  RIVER  IN  NEBRASKA 


This  study  was  motivated  by  a  concern  for  the  habitat  associated  with  the  Platte  River  in 
Nebraska  during  the  winter.  Habitat  studies  of  this  river  have  been  ongoing  for  several  years 
and  have  included  selection  of  representative  study  sites,  morphologic  descriptions -..nd  other 
efforts  to  evaluate  the  present  habitat  and  expected  effects  on  that  habitat  for  various  flows. 
This  study  is  preliminary  because  it  relies  on  existing  data  and  theory  rather  than  on-site 
measurements  during  the  winter.  The  study  is  also  limited  to  an  assessment  of  the  distribution 
of  depths  and  velocities  during  periods  of  ice  cover.  It  is  expected  that  these  results  will 
provide  the  basis  for  data  input  into  simulations  of  the  physical  habitat.  An  example  of  these 
simulations  is  the  Physical  Habitat  Simulation  System  (Milhous  et  al.  1984). 

The  area  of  concern  on  the  Platte  River  extends  from  the  confluence  of  the  North  Platte 
and  the  South  Platte  near  North  Platte.  Nebraska,  to  near  Grand  Island.  Nebraska.  The 
morphology  of  this  reach,  and  subreaches,  is  more  extensively  described  elsewhere  (Bureau 
of  Reclamation  1988).  The  Platte  River  over  this  reach  has  winter  flows  typically  about 
1 000-2000  cfs.  a  slope  of  about  0.00 1 2  ( Fig.  4)  and  a  sand  bed.  The  river  typically  has  a  broad 
dominant  channel,  often  with  islands  and  often  with  smaller  side  channels.  Depths  for  flows 
of  1000-2000  cfs  are  about  1-2  ft  during  summer,  and  mean  velocities  range  from  about  I 
to  2.5  fps.  The  dominant  channel  is  typically  300-600  ft  wide. 

The  scope  of  the  present  study  did  not  include  actual  measurements  during  ice  conditions 
of  the  study  sites  in  1 989.  However.  I  observed  the  ri  veron  22  January  1 989  from  North  Platte 
to  Grand  Island  by  driving  parallel  to  it  and  viewing  it  at  available  bridges.  I  also  viewed  a 
videotape  taken  from  an  overflight  on  12  February  1988.  The  USGS  provided  copies  of 
stream  gauging  notes  for  water  years  1 979. 1 982  and  1 983  for  the  gauge  locations  at  Overton. 
Odessa  and  Grand  Island  (Fig.  4).  These  years  represent  long  and  short  ice  seasons  and  high 
and  low  flows.  While  not  strictly  representativeof  the  morphology  of  the  river(since  gauging 
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sites  are  selected  for  ease  in  discharge  inventory ).  the  data  do  provide  a  basis  for  comparison 
with  the  theory  presented  earlier.  In  particular  it  enables  comparison  of  mean  depths,  mean 
velocities  and  stage  heights  as  a  function  of  flow  discharge  under  both  ice  and  non-ice 
conditions. 

Figure  5  shows  aerial  photographs  of  the  gauging  locations  at  Overton.  Odessa  and  Grand 
Island.  In  the  analysis  that  follows,  data  from  measurements  with  no  ice  or  floating  slush  only 
are  compared  with  data  from  periods  of  icx- conditions.  Data  for  periods  of  ice  conditions  were 
limited  to  cases  where  the  river  was  at  least  70%  ice  covered.  The  effects  of  partial  coverage 
are  discussed  in  Appendix  A. 

(iauge  height  or  stage 

The  most  easily  measured  and  observed  characteristic  of  rivers  is  the  stage,  or  the  level 
of  the  top  surface  of  the  water.  During  periods  of  ice  cover  the  stage  is  higher  than  in  summer 
for  the  same  discharge  due  to  the  flow  resistance  associated  with  the  ice  cover  and  the 
displacement  due  to  the  ice  cover  thickness.  Figure  6  shows  the  observed  gauge  heights  as 
a  function  of  discharge  for  ice  and  non-ice  conditions  for  the  Overton.  Odessa  and  Grand 
Island  gauging  locations.  The  stage-discharge  relationship  is  well  defined  for  non-ice 


Figure  6.  Stage-discharge  observations  with  and  without  an  ice  cover. 
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Figure  6  (cant'd).  Stage-discharge  observations  with  and 
without  an  ice  cover. 


conditions  but  chaotic  for  ice  conditions,  and  it  is  much  higher  for  the  latter.  From  a  habitat 
standpoint  these  high  winter  stages  may  be  significant,  especially  along  the  banks.  However, 
the  irregularity  of  stage-discharge  data  during  winter  obscures  a  regularity  in  flow  depths 
below  the  ice,  described  next. 

Flow  depths 

Figure  7  shows  the  observed  average  flow  depths  as  a  function  of  discharge  for  the 
Overton,  Odessa,  and  Grand  Island  sites  for  both  ice  and  non-ice  conditions.  The  depth  for 
ice  conditions  is  the  average  depth  beneath  the  ice  covers.  The  depths  for  ice  conditions  are 
somewhat  higher  for  the  same  discharge  than  for  non-ice  conditions,  and  there  is  considerable 
scatter  of  the  data  for  both  ice  and  non-ice  conditions.  As  a  rough  estimate,  the  under-ice  flow 
depths  tend  to  be  10-20 %  greaterthan  for  non-ice  conditions  forequivalent  discharges.  From 
Figure  2  this  corresponds  to  a  ratio  of  /(.//^ranging  from  about  0.5  to  0.7.  which,  if  nb-  0.03. 
corresponds  to  an  ice  roughness  of  about  0.02.  This  is  about  the  value  in  Nezhikovskiy 's  data 
(Table  2)  for  an  ice  cover  0.5  m  thick  composed  initially  of  dense  slush. 

Flow  area 

Figure  8  shows  the  cross-sectional  flow  areas  for  ice  and  non-ice  conditions,  again  for  the 
Overton,  Odessa  andGrand  island  sites.  Since  widths  changed  little  in  this  range  of  discharge, 
these  plots  show  the  same  results  as  those  for  depth.  In  short  the  flow  areas  increased  1 0-20% . 

Mean  velocities 

Figure  9  shows  the  mean  velocities  as  a  function  of  discharge  for  ice  and  non-ice 
conditions,  again  for  the  Overton,  Odessa  and  Grand  Island  sites.  If  the  depths  (and  flow 
areas)  increase,  then  the  velocities  decrease,  and  of  course  this  is  what  isobserved.  A  10-20% 
decrease  in  velocity  is  estimated  in  association  with  ice  conditions. 
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Figure  7.  Flaw  depth-discharge  observations  with  and  without  an 
ice  cover. 


Figure  9.  Mean  veloci  ty-ilisi  I  large  observations  with  anil  with¬ 


out  an  we  cover. 


HABITAT  SIMULATIONS 


The  implications  of  these  changes  for  habitat  simulations  for  periods  of  ice  cover  on  the 
Platte  River  reach  considered  are  straightforward.  During  periods  of  nearly  complete  ice 
cover,  the  flow  depths  and  velocities  can  be  calculated  as  for  non-ice  conditions,  and  the 
depths  then  increase  by  10-20%  and  the  velocities  decrease  by  10-20%.  This  corresponds 
approximately  to  a  ratio  of /;.///b=  0.5  to0.7  inetj  lOand  1 1 .  This  is  a  particularly  simple  result, 
which  is  pleasing. 

There  are  some  cautions  that  need  to  be  raised.  First,  the  distributions  of  velocity  across 
the  cross  sections  as  a  function  of  the  distribution  of  depths  has  not  been  examined  in  detail. 
I  expect  they  will  not  be  much  different  for  ice  and  non-ice  conditions.  Second,  all  the  data 
used  here  were  from  gauging  locations  established  by  the  USGS,  and  they  are  at  locations 
with  channel  confinement  by  bridges  and  with  single-channel  river  geometries.  Where  side 
channels  exist,  significant  redistributions  of  flow  may  occur,  particularly  when  one  of  the 
channels  is  ice  free  while  the  other  is  ice  covered.  This  situation  is  analyzed  in  Appendix  B 
for  an  idealized  case  of  two  channels. 

Finally  the  analysis  presented  here  presumes  knowledge  of  whether  or  not  the  flow  is  ice 
covered.  In  turn,  this  depends  on  air  temperatures  and.  in  the  case  of  the  Platte  River,  on 
thermal  inputs  and  magnitudes  of  flows  in  the  river.  These  are  significant,  and  thermal  effects 
were  clear  in  the  data.  Forexample,  while  the  data  for  all  three  sites  covered  the  same  periods, 
there  were  significantly  fewer  observations  of  ice  cover  at  Overton  than  at  Odessa  and  Grand 
Island.  The  video  coverage  of  February  1988  clearly  showed  similar  effects,  with  ice 
conditions  increasing  in  extent  in  the  downstream  direction,  or  conversely  thermal  effects 
increasing  as  one  goes  upstream. 


Application  of  findings  to  PHABSIM 

The  findings  of  this  study  are  centered  on  the  Platte  River  but  procedures  may  be  applied 
to  other  rivers.  The  procedure  is  as  follows. 

Using  the  hydraulic  analyses  contained  in  PHABSIM  (Milhous  et  al.  1 984),  calculate  the 
velocities  Vq  and  depths  Dq  for  the  cross  section  under  consideration  for  non-ice  periods. 
Then,  for  the  periods  when  an  ice  cover  is  present,  the  velocities  l'.  and  depths  D.  associated 
with  ice  conditions  are  given  by  eq  1 1  and  10.  respectively,  repeated  here: 


V-,  =  0.76  P0 


.1  +  (/ii/«b)3/2 


2/5 


(ID 


D,  _ 


Dn 


=  1.32 


1  +  («i/"b)3/2 


12/5 


(10) 


To  use  these  equations,  the  Manning’s  roughness  coefficients  for  the  bed  nband  for  the  ice 
/(.must  be  specified.  Since  »b  must  be  specified  for  the  open-channel  hydraulic  analyses,  it 
is  presumably  already  available.  The  value  of  n.  should  be  specified  by  using  Table  2  as  an 
approximate  guideline  or  by  direct  measurements  in  the  field.  Alternatively  the  combined 
effect  of  n./nb  can  be  determined  directly  by  comparing  velocity  and  depth  ratios  as  a  function 
of  discharge  using  detailed  cross-section  measurements  during  both  open-water  and  ice- 
covered  conditions,  as  was  done  in  this  report  and  further  summarized  in  the  next  section. 


Procedures  for  determining  «./nb 

Using  values  of  n.  and  «b  from  other  studies  is  clearly  less  preferable  than  directly 
determining  values  at  the  site  of  interest.  The  value  of  h  is  ordinarily  determined  from  open- 
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channel  measurements  of  discharge,  hydraulic  radius  (or  depth),  flow  area  and  slope  using 
the  Manning’s  equation  (eq  2).  It  is  desirable  to  obtain  t/^from  periods  when  the  water  is  cold, 
since  temperature  can  affect  the  bed  roughness.  It  is  also  assumed  that  /^stays  constant  when 
the  ice  cover  is  added  although  this  assumption  has  not  been  adequately  validated.  With  »b 
specified  and  with  measurements  during  ice  cover  periods  of  discharge,  flow  areas,  depths 
and  slopes,  the  composite  roughncs.,  «  can  be  determined  from  eq  4  and  n.  can  then  be 
determined  from  eq  5.  The  value  of  /»  does  not  stay  constant  throughout  the  winter.  Newly 
formed  covers  tend  to  have  higher  roughness  values,  which  attenuate  somewhat  during  the 
winter,  only  to  increase  somewhat  just  prior  to  break-up.  The  major  effect  of  the  ice  cover 
is  the  addition  of  the  second  boundary,  and  the  roughness,  while  of  significance  in  extreme 
cases,  is  general  ly  less  influential  than  the  effect  of  this  second  boundary.  For  example,  a  50% 
error  in  the  estimate  of  the  value  of  n.  results  in  only  about  a  20%  error  in  the  estimate  of  depth 
or  velocity. 
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APPENDIX  A:  EFFECT  OF  PARTIAL  ICE  COVERAGE 


INTRODUCTION 

It  is  relatively  common  tor  stationary  ice  to  cover  only  part  of  the  width  of  a  river.  This 
occurs  sometimes  during  formation  when  ice  grows  out  from  the  shore.  It  also  often  occurs 
when  there  is  a  warm  tributary  flow  that  results  in  a  narrow  band  of  open  water  downstream 
of  the  entry  point.  For  these  cases  the  equations  developed  in  the  main  report  may  be  applied 
to  the  open  and  ice-covered  portions  and  the  flows  partitioned  between  the  covered  and  open- 
flow  areas  since  the  energy  slope  is  the  same  for  both.  The  following  analysis  of  partially 
covered  flows  is  essentially  the  same  as  published  by  Calkins  et  al.  ( 1 982),  but  it  also  includes 
an  explicit  result  for  the  limiting  case  of  a  very  thin  ice  cover.  This  explicit  result  allows  the 
effect  of  roughness  of  the  ice  cover  on  the  increased  stage  to  be  evaluated. 


ANALYSIS 


The  analysis  starts  with  the  Manning's  equation: 

l  =  L49_  *2/3  s./2  (A1) 

n 

For  the  open  area,  denoted  by  the  subscript  o , 


ro  =  M9_  D2/3  s  1/2 
"b 

where  »b  =  bed  roughness 
D0  =  flow  depth 
5  =  slope. 

For  the  ice-covered  area,  denoted  by  the  subscript  /. 


(A2) 


v :  =  L49_  3  5  1/2  (A3) 

lie  \  2  / 

where  n  is  the  composite  roughness,  which  depends  on  both  the  bed  roughness  n  and  the 
ice  roughness  and  Dt  is  the  depth  of  flow  under  the  ice.  For  ice  with  a  density  of  p.  in  water 
with  a  density  of  pw  and  for  an  ice  thickness  t  the  stage  H  is 


H  -  Da  -  D\  +  -  t,.  (A4) 

Pv, 

The  coverage  C  has  a  value  between  0  (no  ice)  and  1  (complete  ice  coverage).  The  total 
discharge  Q  is  the  sum  of  the  discharges  in  the  open  and  covered  portions,  so  that,  for  a  total 
width  B, 


Q  =  CBV,D,  +  (1-  C)BV„D„.  <A5) 

Substituting  for  V.  and  l’  results  in 

w  I  O 

q  =  CB^L-  Di  5  m  +  (I  -C)B  D:/'  D,yS  1/2  .  (A6) 

Mi  '  2  /  "b 

Using  eq  A4  this  may  be  expressed  in  terms  of  the  stage  H: 
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Q  =  CB  I**-  - 
»c  ' 


in  -  PL,). 

'  Pw  / 


+  (1  -  C)B  W-V3  5  1/2 
/(b 


For  the  stage  /^associated  with  no  ice  ( C-O ),  eq  A7  becomes 


(2  =  L42_  B//5/3  5  i/2 
"h 


1.49  BSxn- 


Equation  A7,  after  some  algebra,  can  also  be  written  in  the  form 


("  -  ?  f 


P  l _  _  Q" h 

22/3(l -C)nc  1.49  (l-C)BS  1/2 


This  is  the  result  of  Calkins  et  al.  (1982)  and  can  be  solved  for  H  by  iterative  procedures. 

For  the  special  case  of  t  — » 0  (a  rigid  but  infinitely  thin  ice  cover),  eq  A 1 0.  again  after  some 
algebra,  can  be  reduced  to  the  form 

//5/3  _  _ Q»b  »c _  (All) 

1.49  BS  1/2  [22/3  (1  -  C)»c  +  0;b] 

The  stage  relative  to  the  ice-free  conditions,  that  is,  H/Hq.  is 


Wo  (i-C)  +  JZ-ik 

oVi  nc 
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C,  ice  Coverage 


Figure  A I .  Effect  of  partial  cover¬ 
age  hv  a  stationary  ice  cover  on  the 
flow  depth  as  a  function  of  the 
roughness  ratio. 


for  the  limiting  case  of  a  very  thin  ice  cover.  Equation  A 1 2  can  be  put  in  a  form  containing 
/i  //i  explicitly  using  the  Sabaneev  formula,  which  is 


«c 
«  h 


2/3 


( A 1 3 ) 


This  results  in  an  expression  for  H/H  explicitly  in  tenns  of  /».,  n^and  coverage  C  in  the  form 


/  1  .  \ 

1 

+ 

r  i 

2/3 

, 

1  + 1 

b;h! 

?/2 

(A  14) 


Figure  A 1  shows  eq  A 14  for  various  values  of  /i.//ibas  a  function  of  C.  The  value  of  H/H  is 
the  "hydraulic  effect”  of  the  top  boundary  and  roughness.  Actual  H/H  values  should  be 
larger  due  to  the  displacement  effects  of  finite-thickness  ice  covers.  The  case  //.///  =  0  is.  of 
course,  physically  impossible  since  it  represents  an  ice  cover  with  no  roughness  or  friction. 

The  usefulness  of  Figure  A 1  is  that  it  shows  that  the  stage  rise  associated  with  imposition 
of  an  ice  cover  does  not  abruptly  change  as  final  closure  occurs,  but  rather  increases  steadily 
as  the  coverage  goes  from  0  to  1 . 
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APPENDIX  B:  ICE  EFFECTS  ON 
MULTIPLE-CHANNEL  FLOW  DISTRIBUTION 

INTRODUCTION 

It  is  not  uncommon  for  rivers  to  have  multiple  channels,  and  the  distribution  of  the  total 
flow  to  the  individual  channels  depends  on  their  relative  hydraulic  characteristics  of  depth, 
width  and  surface  roughness.  During  winter  the  formation  of  ice  covers  on  these  channels 
may  change  the  distribution  of  the  total  flow.  Here  I  examine  the  extent  of  some  of  these  ice 
effects  on  multiple  channel  flow  distribution  for  the  case  of  a  two-channel  situation. 


ANALYSIS 

1  start  by  imposing  the  end  conditions,  and  I  denote  the  channels  with  subscripts  I  and  2. 
and  characterize  their  widths  B (  and  B,.  their  mean  depths  D  { and  D,,  and  their  discharges 
(?  and  Q..  The  total  discharge  is  Q  and  the  total  width  is  B.  1  assume  that  the  channels  emerge 
from  a  single  channel  and  join  again  at  distance  L  downstream,  with  both  channels  having 
the  same  bed  slope  S  and  equal  open-channel  roughnesses  characterized  by  a  Manning's  nb 
value,  where  the  b  subscript  denotes  the  roughness  associated  with  the  bottom  under  open- 
surface  conditions. 


Figure  Bl  Definition  .sketch  for  the  analysis  of  ice  effects  on  the 
flow  distribution  in  multiple  channels. 


Figure  B I  shows  a  schematic  of  the  case  where  one  channel  is  ice  covered  while  the  other 
is  ice  free.  In  channel  I  the  depth  is  the  distance  from  the  bottom  to  the  water  surface,  while 
in  channel  2  the  depth  of  flow  is 

£>->  =  D\  -  -B-  r  (Bl) 

Pw 

where  p  and  pu  are  the  densities  of  ice  and  water,  and  t  is  the  ice  thickness.  The  velocity  in 
channel  I  is 


\  ,  =  L42_  Df/3  5  1/2  (B2) 

"b 

md  in  the  ice-covered  channel  2. 

IS  =  LW  telP  s  1/2  (B3) 

"c  '  2  ’ 
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Figure  B2.  Relationship  between 
ice-bottom  roughness  ratio  and 
composite-bottom  roughness  ra¬ 
tio  using  the  Sabaneev  formula. 


where  « ,is  the  composite  roughness  of  the  ice-covered  channel.  The  Sabaneev  formula  for 
composite  roughness  gives  n  in  terms  of  the  open-channel  (bottom)  roughness  »hand  the  ice- 
cover  roughness  n: 


(B4) 


As  ranges  from  0.5  to  2,  n Jn^  ranges  correspondingly  from  0.77  to  1.54  (Fig.  B2). 
The  total  discharge  is  composed  of  the  discharges  in  the  two  channels: 


Q  =  Q  i  +  Qj  =  B\  ^  i  D\  +  $2  l7 2  7^2 


=  B,L42-Di5/3<>,/2  +  BiLtH.  OlL  .  (B5) 

nh  ~  nc  22/3 


The  fraction  of  the  discharge  carried  in  channel  1  is 


B\  dP  S  1/2 

Q±  =  _ lb _  (B6) 

Q  IADi73  S  1/2  +  Bt  s  1/2 

«b  '  «c  22/3 


«c  =  «b 


I  + 


(S) 


,p° 


2/3 


which  can  be  written  in  the  form 

Ql  =  _ ! _ .  (B7) 

Q  1  +  g2  «b  ( ,  _  Pi  >  f 
22/3  B]  «c  \  pw  £>i  / 

Clearly,  if  both  channels  are  ice  free,  the  flow  is  distributed  proportional  to  the  widths,  that 
is. 
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The  limiting  case  of  a  thin  ice  cover  corresponds  to  a  skim  of  ice  on  channel  2  such  that 
t/D  — >  0  and 


Q±  =  1 _  .  (B9) 

Q  1  +  ft.  Uh.  _]_ 

fi,  «c  22/3 

This  relationship  is  shown  in  Figure  B3  as  a  function  of  the  alternate  terms  BJB  and  »./nb. 
Figure  B4  shows  the  discharge  carried  by  the  ice-covered  channel  2  relative  to  the  discharge 
carried  by  the  same  channel  with  no  ice  is  shown  (in  both  cases  with  the  larger  channel  ice 
free).  This  ratio  is  denoted  by  Q*. 


Figure  B3.  Effect  of  flow 
in  channel  1  of  a  skim  of 
ice  on  channel  2  as  a  func¬ 
tion  of  roughness  ratio. 


In  Figure  B3  the  effect  of  a  thin  skim  of  ice  on  one  channel  is  to  increase  the  proportion 
of  flow  in  the  ice-free  channel  relative  to  the  “both  channels  ice  free”  case  (equivalent  to  the 
n./«b  =  0  line  in  the  figure).  The  effect  is  not  much  for  the  larger  channels,  ice  free  or  ice 
covered.  In  Figure  B4  the  effect  on  discharge  of  the  smaller  channel  is  shown.  Here  the  effect 
is  significant,  especially  for  a  small  side  channel  covered  with  ice.  As  an  example  for  a  side 
channel  with  a  width  that  is  20%  of  the  total  width  of  both  channels  and  with  a  ratio  of 
M./«b=  1,  the  discharge  in  the  smaller  channel  is  68%  of  what  it  would  be  under  non-ice 
conditions;  adding  finite  thickness  to  the  ice  cover  will  reduce  the  side-channel  discharge 
even  further. 

Similar  results  hold  in  the  opposite  case,  that  is,  a  small  channel  maintained  ice  free  while 
the  wider  main  channel  is  ice  covered.  In  this  case  flow  is  diverted  to  the  smaller  channel. 


Bg/B 


Figure  B4.  Ratio  of  flow  in 
ice-covered  channel  2  (when 
channel  I  is  ice  free)  to  com¬ 
pletely  ice-free  flow  as  a  func¬ 
tion  of  roughness  ratio. 


Figure  B5.  Ratio  of  flow 
in  ice-covered  channel  I 
(when  channel  2  is  ice 
free)  to  completely  ice- 
free  flow  as  a  fur  tion  of 
roughness  ratio. 


These  results  are  shown  in  Figure  B5.  In  this  case  covering  the  main  channel  with  ice  while 
maintaining  the  smaller  channel  ice  free  results  in  a  significant  increase  in  the  discharge  in 
the  smaller  channel. 

These  results  are  not  meant  to  be  definitive  for  actual  field  cases,  which  should  be  analyzed 
using  conventional  hydraulic  analyses  with  ice-cover  effects.  In  some  cases  the  distribution 
of  discharge  is  determined  not  by  ice-cover  effects  as  analyzed  here  but  by  other  hydraulic 
controls.  However,  it  is  clear  that  ice  can  have  significant  effects  on  the  flow  distribution 
among  channels. 
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